Regular coordinate systems for Schwarzschild and other spherical spacetimes 
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The continuation of tiie Schwarzschild metric across the event horizon is a well understood problem 
discussed in most textbooks on general relativity. Among the most popular coordinate systems that 
are regular at the horizon are the Kruskal-Szekeres and Eddington-Finkelstein coordinates. Our 
first objective in this paper is to popularize another set of coordinates, the Painleve-GuUstrand 
coordinates. These were first introduced in the 1920's, and have been periodically rediscovered 
since; they are especially attractive and pedagogically powerful. Our second objective is to provide 
generalizations of these coordinates, first within the specific context of Schwarzschild spacetime, and 
then in the context of more general spherical spacetimes. 
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I. INTRODUCTION 

The difficulties of the Schwarzschild coordinates 
{t,r,9,(t)) at the event horizon of a nonrotating black hole 
provide a vivid illustration of the fact that in general rela- 
tivity, the meaning of the coordinates is not independent 
of the metric tensor ga/s- The Schwarzschild spacetime, 
whose metric is given by (we use geometrized units, so 
that c = G = 1) 



-fd^ + f-Ur^ 



(1.1) 



/ = 1 - 2M/r, 



where 



d02 



sin 6 dcf)'^ , indeed gives one of the 



simplest example of the failure of coordinates which have 
an obvious interpretation in one region of the spacetime 
(the region for which r ^ 2M), but not in another (the 
region for which r < 2M). Understanding this failure 
of the "standard" coordinate system is one of the most 
interesting challenges in the study of general relativity. 
Overcoming this obstacle is one of the most rewarding 
experiences associated with learning the theory. 

Most textbooks on general relativity discuss the 
continuation of the Schwarzschild solution across the 
event horizon either via the Kruskal-Szekeres (KS) co- 
ordinates, or via the Eddington-Finkelstein (EF) coor- 
dinates; both coordinate systems produce a metric that 
is manifestly regular at r = 2M. The main purpose 
of this paper is to show that useful alternatives exist. 
One of them, the Painleve-GuUstand (PG) coordinates, 
are especially simple and attractive, and we will consider 
them in detail. We will also generalize them into a one- 
parameter family of coordinate systems, and show that 
the EF and PG coordinates are members of this family. 

In a pedagogical context, the KS coordinates come 
with several drawbacks. First, the explicit construction 
of the KS coordinates is relatively complicated, and must 
be carried out in a fairly long series of steps. Second, 
the fact that r is only implicitly defined in terms of the 
KS coordinates makes working with them rather difficult. 
Third, the manifold covered by the KS coordinates, with 
its two copies of each surface r = constant, is unnec- 
essarily large for most practical applications; while the 



extension across the event horizon is desirable, the pres- 
ence of another asymptotic region (for which r ^ 2M) 
often is not. While the KS coordinates are not to be dis- 
missed out of hand — they do play an irreplaceable role 
in black-hole physics, and they should never be left out 
of a solid education in general relativity — we would ad- 
vocate, for pedagogical purposes and as a first approach 
to this topic, the construction of simpler coordinate sys- 
tems for extending the Schwarzschild spacetime across 
the event horizon. 

A useful alternative are the EF coordinates {v,r,9,(j)), 
in which the metric takes the form 



ds^ 



-fdv^ + 2dvdr + r^ dW 



(1.2) 



The new time coordinate v is constant on ingoing, radial, 
null geodesies (r decreases, 9 and (j) are constant); it is 
related to the Schwarzschild time thy v = t + r* ^ where 



dr 
1 



2M In 



r 
2M 



1 



(1.3) 



The metric of Eq. (1.2) is regular across the event hori- 
zon. While its nondiagonal structure makes i t sl ightly 
harder to work with than the metric of Eq. (1.1), the 
fact that r appears explicitly as one of the coordinates 
makes it much more convenient than the KS version of 
the Schwarzschild metric. We believe that in a peda- 
gogical context, the Eddington-Finkelstein coordinates 
should be introduced before the KS coordinates. 

Our first objective in this article is to popularize an- 
other set of coordinates for Schwarzschild spacetime, and 
propose this system as a useful alternative to the EF co- 
ordinates. These are the Painleve-Gullstrand (PG |,|) 
coordinates {T,r,9,(j)). They are constructed and dis- 
cussed in Sec. IL Our second objective is to provide gen- 
eralizations of this coordinate system. In Sec. Ill we dis- 
cuss a one-parameter family of PG-like coordinates for 
Schwarzschild spacetime. To the best of our knowledge 
this family was first discovered by Kayll Lake in 1994 
1^, but a related family of coordinates was previously 
discussed by Gautreau and Hoffmann 0-^. We show in 
Sec. Ill that the PG and EF coordinates are both mem- 
bers of Lake's family. In Sec. IV we generalize this fam- 
ily of coordinate systems to other spherical (and static) 
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spacetimes; equivalent coordinates were constructed, in 
a two-dimensional context, by Corley and Jacobson [[lO[ . 
In Sec. V we look back at our coordinates, and offer some 
additional comments regarding their construction. In the 
Appendix we relate these coordinate systems to the KS 
coordinates, and provide details regarding the spacetime 
diagrams of Figs. 1 and 2. 



II. PAINLEVE-GULLSTRAND COORDINATES 

It is often a good strategy, when looking for regular 
coordinate systems, to anchor the coordinates to a spe- 
cific family of freely moving observers . We shall em- 
ploy this strategy throughout this paper. The following 
derivation of the PG coordinates can be found in the book 
by Robertson and Noonan . Other derivations can be 
found in Refs. in which the PG coordinates were 

independently rediscovered. 

We consider observers which move along ingoing, ra- 
dial, timelike geodesies of the Schwarzschild spacetime (r 
decreases, 9 and are constant). It is eas y to check that 
in the standard coordinates of Eq. (LI), the geodesic 
equations can be expressed in first-order form as 



E 
7' 



r2 +f^E^ 



(2.1) 



where an overdot denotes differentiation with respect to 
the observer's proper time, and E — E/m is the ob- 
server's (conserved) energy per unit rest mass. (For a 
derivation, see Chap. 11 of Ref. §, Chap. 25 of Ref. §, 
or Chap. 6 of Ref. |^.) We assume that r < 0, and 
the energy parameter is related to the observer's initial 
velocity Voo — the velocity at r = cx) — by 



E 



(2.2) 



In this section we specialize to the particular family of 
observers characterized hy E = 1; our observers are all 
starting at infinity with a zero initial velocity: Voo = 0. 
For these observers, the geodesic equations reduce to 
t = 1/f and r = — yT^^. We notice that Ua, the co- 
variant components of the observer's four-velocity, whose 
contravariant components are = (t, r, 0, 0), is given by 
Ua — (— 1, — — ///, 0, 0). This means that Ua is equal 
to the gradient of some time function T: 



where 



T = t+ I ^^dr. 



(2.3) 



(2.4) 



Integration of the second term is elementary, and we ob- 
tain 



T = t + 4M ( yJrjlM + - In 



yJrjlM - 1 



^JrjlM + 1 



(2.5) 



This shall be our new time coordinate, and (T, r, 0, 0) 
are nothing but the PG coordinates. It should be clear 
that the key to the construction of the PG coordinates 
is th e fa ct that the four-velocity can be expressed as in 
Eq. (2.3). In Sec. V we will explain how this equation 
comes about. 

Going back to Eq. (2.4), we see that dt — dT — 
f^^ \J2M /r dr. Substituting this into Eq. (1.1) gives 

ds^ = -fdT^ + 2^2M/rdTdr + dr^ + dfl^ . (2.6) 

This is the Schwarzschild metric in the PG coordinates. 
An equivalent way of expressing this is 

ds^ = -dT^ + (dr + yj2M/rdT^^ + dn^. (2.7) 

This metric is manifestly regular at r = 2Af , in corre- 
spondence with the fact that our observers do not con- 
sider this surface to be in any way special. (The metric 
is of course still singular at r — 0.) While the metric 
is now nondiagonal, it has a remarkably simple form. It 
is much simpler than the Kruskal-Szekeres metric, and 
we believe that it provides a useful alternative to the 
Eddington-Finkelstein form of the metric, Eq. 1.2. 

In Fig. 1 we show several surfaces T = constant in 
a Kruskal diagram. The construction is detailed in the 
Appendix. The diagram makes it clear that the PG co- 
ordinates do not extend inside the past horizon of the 
Schwarzschild spacetime — the "white-hole region" is not 
covered. The reason for this is that our observers fall in- 
ward from infinity and end up crossing the future, but 
not the past, horizon. By reversing the motion (choos- 
ing the opposite sign for r), we would obtain alternative 
coordinates that extend within the past horizon but do 
not cover the black-hole region of the spacetime. While 
the PG coordinates do not cover the entire KS mani- 
fold, they do cover the two most interesting regions of 
the maximally extended Schwarzschild spacetime. 

Equations (2.6) and (2.7) reveal the striking property 
that the surfaces T = constant are intrinsically flat: Set- 
ting dT = returns ds-^ = dr^ + dfi^ , which is the 
metric of flat, three-dimensional space in spherical polar 
coordinates. The information about the spacetime cur- 
vature is therefore entirely encoded in the "shift vector" , 
the off-diagonal component of the metric tensor. We con- 
sider this aspect of the PG coordinates to be their most 
interesting property. 

We note that it is possible to construct PG-like coor- 
dinates for the nonspherical Kerr spacetime. This was 
carried out by C. Doran in a recent paper |14|. 



III. GENERALIZATION TO OTHER OBSERVERS 

It is easy to generalize the preceding construction to 
other families of freely moving observers. In this sec- 
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FIG. 1. Surfaces T = constant in a Kruskal diagram. The 
vertical and horizontal axes correspond to the Kruskal r and p 
coordinates, respectively; these are defined in the Appendix. 
The thick, diagonal lines represent the two copies of the sur- 
face r = 2M; the future horizon is oriented at +45 degrees, 
while the past horizon is oriented at —45 degrees. The thick, 
hyperbolic line represents the curvature singularity at r = 0. 
The thin lines represent the surfaces T = constant. From the 
bottom up we display the surfaces T = -2M, T = 0, T = 2M, 
T = 3M, T = AM, T = 5Af, T = 6M, and T = 7M. 



tion we consider families such that E is the same for all 
observers within the family, but not equal to unity (as 
in Sec. II). Each family is therefore characterized by its 
unique value of the energy parameter. We find it con- 
venient to use instead the parameter p, related to the 
energy and initial-velocity parameters by 



(3.1) 



We take p to be in the interval < p < 1, with p — 1 
taking us back to the PG coordinates To each value 
of p in this interval corresponds a family of observers, 
and a distinct coordinate system. We are therefore con- 
structing a one-parameter family of PG-like coordinates 
for Schwarzschild spacetime. 

With the geodesic equations now given by i = l/(y^/) 
and f — —y/l — pf / y/p, we find that Ua is now equal to 
a constant times the gradient of a time function T: 



with 



T = t 



f 



dr. 



(3.2) 



(3.3) 



Integration of the second term doesn't present any essen- 
tial difficulties, and we obtain 



T = t-)-2A/^i-^+ln 



FIG. 2. Surfaces T — constant in a Kruskal diagram. The 
axes and the meaning of the thick lines are explained in the 
caption of Fig. 1. The eight bundles of thin lines refer to the 
same values of T as in Fig. 1, from T — —2M (bottom bundle) 
to T = 7M (top bundle). Within a single bundle, each of 
the four lines come with a distinct value of the parameter 
p. The solid line corresponds to p = 1, the long-dashed line 
corresponds to p = 3/4, the short-dashed line corresponds to 
p = 1/2, and the dotted line corresponds to p = 1/4. 



■p/2 



In 



(3.4) 



This shall be our new time coordinate. In Sec. V we will 
return to the question of the origin of Eq. ( p^ ) . 

With dt now equal to dT~ f~^y/l — pf dr, we find that 
the Schwarzschild metric takes the form 



ds^ = -f dT^ + 2^1- pf dTdr +pdr^ + 



(3.5) 



ds' 



1 



1 



dT^ +p{ dr + -Jl -pfdT 
P \ P 



(3.6) 



This metric is still regular at r = 2M, although it is now 
slightly more complicated than the PG form. 

In Fig. 2 we show several surfaces T = constant in a 
Kruskal diagram, for several values of p. This construc- 
tion is detailed in the Appendix. 

In this generalization of the PG coordinates, the sur- 
faces T — constant are no longer intrinsically flat. The 
induced metric is now ds^ = pdr^ + dVl^ , and al- 
though the factor of p in front of dr^ looks innocuous, 
it is sufficient to produce a curvature. It may indeed 
be checked that the Riemann tensor associated with this 
metric is nonzero. The only nonvanising component is 
R% = -(1 -P)/P: and R-'-^Rabcd = A{1 ~ pf / {pr^ f . 

It is instructive to go back to Eq. ( |3.4| ) and check 
that in the limit p — *■ 1, T reduces to the expression 
of Eq. ( |2.5| ). (This must be done as a limiting procedure, 
because T is ambiguous for p = 1.) Taking the limit gives 
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limT = t + 2M , + In 



(3.7) 



which is indeed equivalent to Eq. ( |2.5D . The PG coordi- 
nates are therefore a member of our one-parameter fam- 

iiy- 

Another interesting Hmit is p — > 0, which corresponds 
to — > oo, or Voo — > 1. In this hmit, our observers start 
at infinity with a velocity nearly equal to the speed of 
light. Starting from Eq. (3.4) we have 



In 



/ 



l-f 



t + r* 



(3.8) 



where we have compared with Eq. (1.3). Thus, T = v 
in the limit p 0, and our generalized PG coordi- 
nates reduce to the Eddington-Finkelstein coordinates of 
Eq. ( |1.2| ). This is not entirely surprising, in view of the 
fact that our observers become light-like in this limit. 
The EE coordinates, therefore, are also a (limiting) mem- 
ber of our one-parameter family. 

We have constructed an interpolating family of coordi- 
nate systems for Schwarzschild spacetime; as the parame- 
ter p varies from 1 to 0, the coordinates go smoothly from 
the Painleve-GuUstrand coordinates to the Eddington- 
Finkelstein coordinates. This one-parameter family of 
coordinate systems was first discovered by Kayll Lake 
1^, but a related family of coordinates, corresponding 
to p > 1, were previously introduced by Gautreau and 
Hoffmann Lake obtained the new coordinates by 

solving the Einstein field equations for a vacuum, spher- 
ical spacetime in a coordinate system involving r and 
an arbitrary time T. The intimate relation between his 
coordinates and our families of freely moving observers 
remained unnoticed by him. 



IV. GENERALIZATION TO OTHER 
SPACETIMES 

The coordinates constructed in the previous two sec- 
tions can be generalized to other static and spherically 
symmetric spacetimes. In the usual Schwarzschild coor- 
dinates, we write the metric as 



(4.1) 



where / and '0 are two arbitrary functions of r. Under 
the stated symmetries, Eq. (4.1) gives the most general 
form for the metric. We assume that the spacetime is 
asymptotically flat, so that / — > 1 and ^ as r — > 
oo. If the spacetime possesses a regular event horizon at 
r = To, then /(ro) = and ■0 must be nonsingular for all 
values of r 7^ 0. 

The geodesic equations are now 



E 



(4.2) 



where E is still the conserved energy per unit rest mass. 
Re-introducting p — 1/E'^, we find that the covariant 
componen ts o f the four-velocity can be again expressed 
as in Eq. ( |3.2|) , with a time function T now given by 



T = t + 



Pf 



f 



dr. 



(4.3) 



The second term can be integ rate d if / and ip are known. 
Rewriting the metric of Eq. (4.1) in terms of dT yields 



ds^ = -fe^''' dT^ ^ 
s^^ dr^ 



' pe 



2e^'^^/e~^'l' -pf dTdr 
r^dif, (4.4) 



ds^ = -- dT^ + pe'^'^ (dr+- ^e-2V' - pf dT 
p \ P 



r2 dn^. 



(4.5) 



This metric is manifestly regular at an eventual event 
horizon, at which / vanishes. 

The surfaces T = constant have an induced metric 
given by ds^ = pe^^ dr"^ + r'^ dVL^ . Unless ^ — Q and 
p = I, these surfaces are not intrinsically flat Q. 



V. FINAL COMMENTS 

In all the cases considered in Sees. II, III, and IV, the 
construction of our coordinate systems relied on the key 
fact that the four-velocity could be expressed as Ua — 
—daT/ with p a constant. [This is Eq. (3^), and 
in Sec. II, p was set equal to unity] This property is 
remarkable, and it seems to follow quite accidentally from 
the equations of motion. There is of course no accident, 
but the point remains that not every four-velocity vector 
can be expressed in this form. 

A standard theorem of differential geometry (for ex- 
ample, see Appendix B of Ref. ||]) states that for Ua to 
admit the form of Eq. ( |3.2| ), it must satisfy the equa- 
tions u'^.pU^ = and M[a;/3U^] = 0, in which a semicolon 
denotes covariant differentiation and the square brack- 
ets indicate complete antisymmetrization of the indices. 
The second equation states that the world lines are ev- 
erywhere orthogonal to a family of spacelike hypersur- 
faces, the surfaces of constant T. This ensures that the 
four- velocity can be expressed as Ua — daT , for some 
function /i(a;"). In general, this fun ction is not a con- 
stant, and we do yet have Eq. (3.2). For this we need 
to impose also the first equation, which states that the 
motion is geodesic. When both equatio ns h old we find 
that /i = constant, and this gives us Eq. (3.2). 

In our constructions, we have enforced the geodesic 
equation by selecting freely moving observers. By se- 



lecting radial observers, we have also enforced the condi- 
tion that the geodesies be hypersurface orthogonal. Our 
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strategy for constructing coordinate systems is there- 
fore limited to radial, freely moving observers in static, 
spherically-symmetric spacetimes; it may not work for 
more general motions and/or more general spacetimes. 



of Figs. 1 and 2 are then produced by switching to the 
coordinates t and f and plotting the parametric curves. 
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APPENDIX A: KRUSKAL DIAGRAMS 

The Kruskal diagrams of Figs. 1 and 2 are constructed 
as follows. 

From the Schwarzschild coordinates t and r we define 
two null coordinates, u = t—r* and v — i-f-r*, where r* is 



given by Eq. (1.3). From these we form the null KS coor- 
dinates, V = e^*^ and U = =Fe~'"/''*^, where the upper 
sign refers to the region r > 2M of the Schwarzschild 
spacetime, while the lower sign refers to r < 2M. From 
this we derive the relations 



r 

\2M 



1 



and 



V 
U 



(Al) 



(A2) 



Timelike and spacelike KS coordinates are then defined 
hy V = T + p and U — t — p. In our spacetime diagrams, 
the T axis runs vertically, while the p axis runs horizon- 
tally. The future horizon is located at J7 = 0, and the 
past horizon is at y = 0. The curvature singularity is 
located at UV = 1. 



We express the time function of Eq. (3.4) as 



T = t 



Sir), 



(A3) 



where S{r) is the function of r that results when the 
seco nd t erm of Eq. ( ^.4[ ) is shifted by — r*, as given in 
Eq. (1.3); this function is regular at r = 2M. With this 
definition we have v^T— S,u~T — S — 2r*, as well as 



and 



V = e^/^M^-s/4M (A4) 



U = -e''/2*f {r/2M - l)e-^/4Mg5/4Af ^ (^5) 



The surfaces T = constant give rise to parametric equa- 
tions of t he for m V jr) and U{r), which are obtained from 
Eqs. (A4) and ( |A5| ) by explicitly evaluating the function 
S(r). In these equations, r can be varied from zero to an 
arbitrarily large value without difficulty. The diagrams 
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